Pairing mechanism of high-temperature superconductivity: Experimental constraints 
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Developing a theory of high-temperature superconductivity in copper oxides is one of the out- 
standing problems in physics. It is a challenge that has defeated theoretical physicists for more 
than twenty years. Attempts to understand this problem are hindered by the subtle interplay 
among a few mechanisms and the presence of several nearly degenerate and competing phases in 
these systems. Here we present some crucial experiments that place essential constraints on the 
pairing mechanism of high-temperature superconductivity. The observed unconventional oxygen- 
isotope effects in cuprates have clearly shown strong electron-phonon interactions and the existence 
of polarons and/or bipolarons. Angle-resolved photoemission and tunneling spectra have provided 
direct evidence for strong coupling to multiple-phonon modes. In contrast, these spectra do not 
show strong coupling features expected for magnetic resonance modes. Angle-resolved photoemis- 
sion spectra and the oxygen-isotope effect on the antiferromagnetic exchange energy J in undoped 
parent compounds consistently show that the polaron binding energy is about 2 eV, which is over 
one order of magnitude larger than J = 0.14 eV. The normal-state spin-susceptibility data of hole- 
doped cuprates indicate that intersite bipolarons are the dominant charge carriers in the underdoped 
region while the component of Fermi-liquid-like polarons is dominant in the overdoped region. All 
the experiments to test the gap or order-parameter symmetry consistently demonstrate that the 
intrinsic gap (pairing) symmetry for the Fermi-liquid-like component is anisotropic s-wave and the 
order-parameter symmetry of the Bose-Einstein condensation of bipolarons is d-wave. 



BCS THEORY AND THE CONVENTIONAL 
ISOTOPE EFFECT ON Tc 

In 1911 curiosity concerning the electrical properties of 
metals at low temperatures led the Dutch physicist, H. 
K. Onnes and his assistant G. Hoist to discover super- 
conductivity at 4.2 K in mercury [1]. This discovery was 
one of the most important experimental findings in low 
temperature physics. Since then tremendous theoretical 
and experimental efforts have been made with the aim 
of clarifying the microscopic mechanism responsible for 
this novel ground state. 

On the long way towards a microscopic understanding 
of superconductivity, the observation of an isotope effect 
on Tc in 1950 0, Q gave important clues to the under- 
standing of the microscopic theory of superconductivity. 
The presence of an isotope effect thus implies that su- 
perconductivity is not of purely electronic origin. In the 
same year H. Frohlich 3 pointed out that the electron- 
phonon interaction gave rise to an indirect attractive in- 
teraction between electrons, which might be responsible 
for superconductivity. Frohlich's theory played a decisive 
role in establishing the correct mechanism. Later, in 1956 
L. Cooper Q discovered that electrons with an attrac- 
tive interaction form bound pairs (so called Cooper pairs) 
which lead to superconductivity. However, the existence 
of electron pairs does not necessarily imply a phonon me- 
diated pairing. Indeed, Bose condensation as considered 
in 1955 by Schafroth ,(|, is also a possible mechanism 
for superconductivity, but the model was not able to ex- 
plain the existence of an isotope effect. Finally, in 1957, 
Bardeen, Cooper and Schrieffer 7] developed the BCS 
theory that was the first successful microscopic theory of 
superconductivity. 



It is remarkable that the BCS theory shows that 
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where w is a typical phonon frequency (e.g., the Debye 
frequency lud)- The electron-phonon coupling constant 
N{0)V is the product of an electron-phonon interaction 
constant V and the electronic density of states at the 
Fermi surface N{0), both of which are assumed to be 
independent of the ion mass M. 

This formula implies an isotope-mass dependence of 
Tc, with an isotope-effect exponent a = —dlnTc/dlnM 
= 1/2, in excellent agreement with the reported isotope 
effects in the non-transition metal superconductors (e.g., 
Hg, Sn and Pb). On the other hand, many supercon- 
ducting transition metals show a much smaller isotope 
effect (e.g. Zr, Ru: a — ± 0.05), and in some materi- 
als even negative isotope shifts have been reported (e.g. 
U: a = —2.2 ± 0.2). The obvious discrepancy between 
the observed isotope effect and the prediction of the BCS 
theory demonstrates the limitations of the simplified BCS 
approach. Nevertheless, the phonon mediated BCS the- 
ory was the basis for further and more sophisticated mi- 
croscopic theories of superconductivity. 

Since the original BCS theory of 1957 there has been 
considerable progress toward a deeper understanding of 
the electron-phonon interaction. In the strong-coupling 
Eliashberg model 's'l the electron-phonon coupling con- 
stant Aep is related to the spectral function a'^(uj)F(uj), as 
defined by the product of the frequency-dependent aver- 
age electron-phonon interaction^ a^(cj) and phonon den- 
sity of states F{lu). McMillan !9] numerically solved the 
Eliashberg equations and found an expression for Tc tak- 
ing into account the Coulomb repulsion between electrons 
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and the retarded nature of the phonon-induced interac- 
tion 
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where /i* is an effective Coulomb repulsion. The BCS re- 
sult of Eq. ([T]) is recovered for Aep ^ 1, and in this weak- 
coupling limit Aep — substitutes the role of N{0)V. As 
a consequence, in strong-coupling superconductors the 
isotope effect exponent a is no longer a universal value 
and given by 



^ ^ 1 A _ 1.04(1 + 0.62Aep)(l + Aep)M*' 



2V' [Aep-M*(l + 0.62Aep)]2 



(3) 



Within the framework of the Eliashberg theory a small 
or even negative isotope exponent is possible. Therefore, 
the observation of a small isotope exponent could still be 
consistent with the phonon-mediated pairing mechanism. 



the occurrence of superconductivity. For a conclusive 
decision on the role of lattice vibrations in superconduc- 
tivity, one should not only study the isotope effect on 
Tc, but also the isotope effects on some other parameters 
such as the effective supercarrier mass m* and the su- 
percarrier density rig. The latter effects may offer more 
information about the role of phonons, and will be crucial 
to the understanding of the physics of high-temperature 
superconductors. 

Studies of isotope shifts of Tc have been carried out 
in almost all known cuprates. A comprehensive review 
is given by Franck ll|. Most of the studies reported 
so far concern the oxygen isotope shift by replacing -'^^O 
with ^^O, partly because the experimental procedures are 
simple and reliable. Today it is generally accepted that 
the cuprates doped for maximum Tc exhibit a small but 
clearly non- vanishing ao- This is true for all the different 
families of doped cuprates. Some of the results arc listed 
in Table. HI 



ISOTOPE EFFECT ON n IN HIGH-Tc CUPRATE 
SUPERCONDUCTORS 

Since the discovery of superconductivity in 
La2_a;Baj;Cu04 by Bednorz and Miiller in 1986 
[Tq] . tremendous efforts have been made to clarify the 
microscopic pairing mechanism for high-temeprature 
superconductivity. In some of the cuprates the critical 
temperatures are far beyond a supposed maximum 
Tc « 30 K, estimated theoretically on the basis of the 
conventional phonon-mediated mechanism fo*]. However, 
the theoretical estimate has not been well justified. 

Contrary to low Tc elemental superconductors, the 
cuprates have a complex crystal structure and consist 
of at least three constituents, leading to complicated 
phononic and electronic densities of states. The struc- 
ture of cuprates may be described by the three main fea- 
tures: i) the Cu02 planes, ii) the apical oxygen, iii) the 
metal-oxygen charge reservoir. All the high-Tc super- 
conductors have Cu02 planes, but the number of planes 
varies among the different families of cuprates, ranging 
from a single-layer up to an infinite-layer structure. 

The investigation of the isotope effect in cuprates is 
certainly much more difficult than in elemental super- 
conductors. In compound superconductors, such as the 
cuprates, it is convenient to define a partial isotope ex- 
ponent for a constituent with mass Mi as 



Tc (K) Oxygen-isotope exponent ao 
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where the index i denotes a constituent with mass Mi 
(e.g. oxygen). The magnitude of ai may provide insight 
into which phonons are important in the occurrence of su- 
perconductivity. As discussed above, small or even van- 
ishing partial isotope exponents do not necessarily imply 
that the corresponding constituent is not important for 



Lai.85Sro.i5Cu04 


35-37 


0.08-0.13 


YBa2CU306.94 


90-92 


0.02-0.03 


YBa2Cu408 


79.5 


0.06-0.08 


Bi2Sr2CaCu2 08+!; 


75 


0.03-0.05 


Tl2Ca2BaiCu30io 


121 


±0.12 


Ndi.ssCeo.isCuOi 


24.5 


< 0.05 



TABLE I: Oxygen isotope effect for various cuprates doped 
for optimum Tc (compiled from Ref. ^llj]) . 



Although the different compounds exhibit quite differ- 
ent Tc's, it is striking that the oxygen-isotope exponents 
(OIEs) are similar in magnitude and remain small. This 
suggests that the OIE depends neither on Tc nor on spe- 
cific structural properties such as the number of Cu02 
planes and the electronic properties of the metal-oxygen 
buffer layers. The small OIE observed in the optimally 
doped cuprates suggest that phonons might not be im- 
portant in high temperature superconductivity. This has 
led to many exotic n on-p honon mediated pairing mecha- 
nisms (e.g., see Ref. [12|). 

However, the doping dependence of the OIE has been 



extensively studied in different cuprate systems isl- 19j| 
For a particular family of doped cuprates the OIE in- 
creases with decreasing Tc, and can be even larger than 
the BCS prediction 13^. Such an anomalously large OIE 
may imply that phonons play an important role in the 
occurrence of high temperature superconductivity in the 
copper oxides. Any correct theories for describing the 
physics of high- Tc superconductors must consistently ex- 
plain both the small OIE in the optimally-doped samples, 
and the anomalously large OIE in the underdoped sam- 
ples. 
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The non-vanishing OIEs indicate that the oxygen- 
dominated phonon modes are involved in the occurrence 
of superconductivity in the cuprates. As we will see 
below, there are strong coupling features to multiple 
phonon modes in tunneling and angle-resolved photoe- 
mission spectra. This implies that there should be iso- 
tope shifts related to other atoms. Indeed, large copper- 
isotope shifts, with acu being close to ao, have been 
observed in underdoped La2_a;Sra;Cu04 (LSCO) [l^ and 
in oxygen-depleted YBa2Cu307_y [2^. This implies that 
Cu-dominated phonon modes also play important role in 
the pairing. 



OXYGEN-ISOTOPE EFFECT ON THE IN-PLANE 
SUPERCARRIER MASS IN CUPRATES 



The conventional phonon-mediated superconducting 
theory is based on the Migdal adiabatic approximation 
in which the phonon-induced electron self-energy is given 
correctly to order (m/M)^/^ ~ 10~^, where m is the 
mass of an electron. Within this approximation, the 
density of states at the Fermi level N(0), the electron- 
phonon coupling constant Aep, and the effective mass of 
the supercarriers are all independent of the ion- mass M . 
This approximation has been working very well in con- 
ventional elemental superconductors because the ratio of 
the sound velocity is much smaller than the Fermi ve- 
locity. However, since the phase velocity of long wave- 
length optical phonons is large compared to the Fermi 
velocity, it may break down in some compounds with 

The 



carriers m* = mexp{'~fEp/hu!), this definition leads to 



/3 = — — ln(m*/m). 



(5) 



reasonably strong electron-phonon interactions [21 



break-down of the Migdal approximation may lead to 
the formation of polarons— new quasi-particles moved to- 
gether with local lattice distortions. The effective mass 
of polarons m* will depend on M. This is because the po- 
laron mass m* = m exjp('jEp/huj) 22, 2^, where m is the 
band mass in the absence of the electron-phonon inter- 
action, 7 is a constant, and w is a characteristic phonon 
frequency which depends on the masses of ions. Hence, 
there should be a large isotope effect on the polaronic 
carrier mass, in contrast to the zero isotope effect on the 
effective carrier mass in ordinary metals. Due to the po- 
laronic band narrow effect, the effective density of states 
at the Fermi level A^(0) and the effective electron-phonon 
coupling constant Aep [proportional to Af(0)] should also 
depend on the ion mass M . This will lead to more com- 
plicated isotope effect on Tc- Therefore, negative, van- 
ishing, or small isotope exponents a ^ 1/2 in optimally 
doped cuprates do not necessarily imply that phonons 
are not important for the pairing mechanism if polarons 
are bound into the Cooper pairs. 

The total exponent of the isotope effect on the effective 
carrier mass is defined as /3 = ^ —dlnm*/d\nMi {Mi is 
the mass of the ith atom in a unit cell). For polaronic 



It is interesting that the simple relation m* = 
m exp{'y Ep / ?iuj) is even valid in the weak coupling region 
in the case of the long-range Frohlich electron-phonon 
interaction [2^. Then the polaron mass enhancement 
factor m*/m in this case is simply equal to exp(— 2/3). 
Therefore, if electron-phonon coupling in a solid is strong 
enough to form polarons and/or bipolarons, one should 
expect a substantial isotope effect on the effective mass 
of carriers. 

An important and essential proof for the existence of 
polarons and/or bipolarons in cuprates is provided by the 
observation of the substantial oxygen-isotope effect on 
the penetration depth [I6, 1^ 24-31|. Zhao et al. made 
the first observation of this effect in optimally doped 
YBa2Cu3 06.93 in 1995 24]. By precisely measuring the 
diamagnetic signals near Tc for the ^^O and ^^O sam- 
ples, the authors were able to deduce the oxygen-isotope 
effect on the penetration depth A(0). It turns out that 
AA(0)/A(0) = 3.2 % (Ref. 0)- 

In fact, for highly anisotropic materials, the observed 
isotope effect on the angle-averaged A(0) is the same 
as the isotope effect on the in-plane penetration depth 
Aab(O). From the magnetic data for YBa2Cu3 06.93, 
Lai.85Sro.i5Cu04, and Bii.6Pbo.4Sr2Ca2Cu30io+j^, one 
finds that AAab(0)/Aa6(0) = 3.2±0.7% for the three op- 
timally doped cuprates ^2£]. Several independent ex- 
periments have consistently shown that the carrier den- 
sities of the two isotope samples are the same within 
0.0004 per unit ceh [11, [23, [s^ ■ Therefore, the ob- 
served oxygen-isotope effect on the in-plane penetration 
depth is caused only by the isotope dependence of the in- 
plane supercarrier mass. Recently, direct measurements 
of the in-plane penetration depth by low energy muon- 
spin-relaxation (LE/iSR) technique [30| have confirmed 
the earlier isotope-effect results. It is found that [s^ 
AAa6(0)/Aa6(0) = 2.8±1.0%. It is remarkable that the 
isotope effect obtained from the most advanced and ex- 
pensive technology (LE/iSR) is the same as that deduced 
from simple magnetic measurements. 

Fig. 1 shows the doping dependence of the exponent 
iPo) of the oxygen- isotope effect on the in-plane super- 
carrier mass in La2-xSra;Cu04. Here the exponent is 
defined as /3o = —dhim'^l/dhiMo- It is apparent that 
the exponent increases with decreasing doping, in agree- 
ment with the fact that doping reduces electron-phonon 
coupling due to screening of charge carriers. The large 
oxygen-isotope effect on the in-plane supercarrier mass 
cannot be explained within the conventional phonon- 
mediated pairing mechanism where the effective mass of 
supercarriers is independent of the isotope mass (3^ . In 
particular, the substantial oxygen-isotope effect on m*^ 
in optimally doped cuprates indicates that the polaronic 
effect is not vanished in the optimal doping regime where 
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FIG. 1: The doping dependence of the exponent Po of 
the oxygen-isotope effect on the in-plane supercarrier mass 
in La2-2!Sr2,Cu04. The exponent is defined as Bp = 
-dlnmll/dlnMo. The data are from Ref. [2i,[2i,[2g| 



the BCS-like superconducting transition occurs. This 
suggests that polaronic carriers may be bound into the 
Cooper pairs in optimally doped and overdoped cuprates. 
This scenario 



22l |29| can naturally explain a small iso- 



tope effect on Tc in optimally doped cuprates. 



have similar elcctron-phonon coupling strengths. Here 
we will show that the small isotope shift of T/v is also 
consistent with a large polaron binding energy due to 
the fact that the electron-phonon interaction renormal- 
izes the antiferromagnetic exchange energy J through the 
forth order process [3^ . 

From the oxygen- isotope shift of Tat, one can deter- 
mine the oxygen-isotope effect on the antiferromagnetic 
exchange energy J, that is, / J ~ —0.6% (Ref. ^i^). 
A slightly larger oxygen-isotope effect on J {A J/ J ~ 
—0.96%) has also been extracted from the mid-infrared 
two-magnon absorption spectra of the oxygen-isotope ex- 
changed YBa2Cu306.o (Ref. [36])- This unconventional 
isotope effect on J can be explained in terms of strong 
electron-phonon coupling and polaronic effect. Eremin et 
al. [3^ have considered strong electron-phonon coupling 
within a three-band Hubbard model appropriate for the 
charge-transfer insulators. They show that the antifer- 
romagnetic exchange energy J depends on the polaron 
binding energy due to oxygen vibrations, on the po- 
laron binding energy E'^^ due to copper vibrations, and 
on their respective vibration frequencies ujq and ucu ■ At 
low temperatures, J is given by [35| 



J= ./o(l 



3E"hLJ, 
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Cu , 



(6) 



DETERMINATION OF THE 
ELECTRON-PHONON COUPLING STRENGTH 
IN UNDOPED PARENT CUPRATES 

The observed substantial oxygen-isotope effect on the 
effective supercarrier mass suggests that doped holes are 
strongly coupled with phonons to form polarons and/or 
bipolarons. It is natural that undoped parent compounds 
like La2Cu04 should have stronger electron-phonon in- 
teractions because no doped charge carriers can screen 
the electron-phonon interactions. Very strong electron- 
phonon interactions can reduce electron hopping integral 
t through the polaronic effect 22], which may lead to 
the renormalization of magnetic exchange energy. Then 
one might expect that the antiferromagnetic (AF) ex- 
change energy in antiferromagnetic cuprates and the fer- 
romagnetic exchange energy in ferromagnetic mangan- 
ites should depend on the isotope mass if both systems 
have very strong electron-phonon coupling. Following 
this simple argument, Zhao and his co-workers initiated 
studies of the oxygen-isotope effect on the AF ordering 
temperature in several parent cuprates 33[ and the 
oxygen-isotope effect on the Curie temperature in ferro- 
magnetic managnites [s^]- A small oxygen- isotope shift 
of Tn was consistently observed in La2Cu04 while a 
giant oxygen-isotope shift of the Curie temperature was 
observed in doped manganites |34l |. The large difference 



in the isotope effects on the magnetic ordering tempera- 
tures appears contradictory since both materials should 



Here Ap^ is the charge-transfer gap, which is measured 
to be about 1.5 eV in undoped cuprates. From the above 
equation, one can clearly see that there should be an ob- 
servable oxygen-isotope effect on J if the polaron binding 
energy is significant. The oxygen-isotope effect on J can 
be readily deduced from Eq. [51 
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The charge-transfer gaps Apd have been measured for 
both La2Cu04 and YBa2Cu306 systems [sll, that is, 
Apd = 1.81 eV for La2Cu04 and Apd = 1.60 eV for 
YBa2Cu306. If we take hujo — 0.075 eV and substitute 
the above parameters into Eq.[7l we obtain = 1.46 eV 
for La2Cu04 and E^ = 1.82 eV for YBa2Cu306. Since 
the copper-related phonon modes have much lower ener- 
gies than the oxygen-related modes, the copper-related 
phonon modes may contribute little to the polaron for- 
mation, that is, i?^" << 1. Therefore, the total polaron 
binding energy Ep should be about 1.5 eV and 1.9 eV for 
La2Cu04 and YBa2Cu306, respectively. 

Angle-resolved photoemission spectroscopy (ARPES) 
data of undoped La2Cu04 have been explained in terms 
of polaronic coupling between phonons and charge car- 
riers [3^ . From the width of the phonon side band in 
the ARPES spectra, the authors find the polaron bind- 
ing energy to be about 1.92 eV, in good agreement with 
their theoretical calculation based on a shell model [SSj. 
On the other hand, the observed binding energy of the 
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side band should be consistent with a polaron binding 
energy of about 1.2 eV (Ref. [si]). This should be the 
lower limit because the binding energy of the side band 
decreases rapidly with doping and because the sample 
may be lightly doped [s!]. Therefore, the ARPES data 
suggest that 1.2 eV < Ep < 1.9 eV, which is in quanti- 
tative agreement with the value {Ep > 1.5 eV) deduced 
from the isotope effect on the exchange energy. 

The parameter-free estimate of the polaron binding en- 
ergy due to the long-range Frohlich-type electron-phonon 
interaction has been made for many oxides including 
cuprates and manganites fs^. The polaron binding en- 
ergy due to the long-range Frohlich-type electron-phonon 
interaction is estimated to be about 0.65 eV for La2Cu04 
(Ref. [3^). The polaron binding energy due to the Qi- 
type Jahn- Teller distortion is about 1.2 eV for La2Cu04 
(Ref. [4o| ) . The total polaron binding energy Ep should 
be at least 1.85 eV, in excellent agreement with the value 
( 1.5 eV < Ep < 1.9 eV) deduced independently from the 
oxygen-isotope effect on J and ARPES data. 

If there are very small amounts of charged carriers in 
these nearly undoped compounds, the optical conduc- 
tivity due to the polaronic charge carriers will show a 
broad peak at Em = ^'jEp (Ref.[39|), where 7 is 0.2—0.3 
for layered cuprates [39|. There appears to exist the 
third broad peak at 0.7-0.8 eV in the optical conductiv- 
ity of SraCuOaCla (Ref. 0), La2Cu04 (Ref. [13]), and 
YBa2Cu306 (Ref. [43|). This peak should be caused by 
the polaronic effect because the energy scale for the peak 
is similar to that predicted from the polaron theory as- 
suming 7 ~ 0.2. Hole doping will reduce the value of Ep 
and thus E,n due to screening of charged carriers. Indeed, 
Em was found to be about 0.6 eV for Lai.98Sro.o2Cu04 
and 0.44 eV for Lai.94Sro.o6Cu04 (Ref. 



STRONG COUPLING TO MULTIPLE-PHONON 
MODES IN ANGLE-RESOLVED 
PHOTOEMISSION SPECTROSCOPY AND 
TUNNELING SPECTRA 

It is known that the electron self-energy renormaliza- 
tion effect in the form of a "kink" in the band dispersion 
can reveal coupling of electrons with phonon modes. The 
"kink" feature at an energy of about 70 meV has been 
seen in the band dispersion of various cuprate supercon- 
ductors along the diagonal ( "nodal" ) direction [45| . The 
energy scale of about 70 meV appears to provide evidence 
for strong coupling between electrons and the 70 meV 
Cu-0 half-breathing mode observed by neutron scatter- 
From the measured dispersion, one can extract 



mg 

the real part of electron self-energy that contains infor- 
mation about coupling of electrons with collective boson 
modes. The fine coupling structures can be revealed in 
the high-resolution ARPES data. Very recently, such 
fine coupling structures have been clearly seen in the 
raw data of electron self-energy of deeply underdoped 



La2_2:Srj;Cu04 along the diagonal direction [47|. Using 
the maximum entropy method (MEM) procedure, they 
are able to extract the electron-phonon spectral density 
a'^F{uj) that contains coupling features at 27 meV, 45 
meV, 61 meV and 75 meV. These ARPES data and exclu- 
sive data analysis 47| clearly indicate that the phonons 
rather than the magnetic collective modes are responsible 
for the electron self-energy effect. 

In fact, the fine coupling structures also appear in the 
earlier high-resolution ARPES data of a slightly over- 
doped BSCCO with Tc = 91 K (Ref. @). Zhao ^ rean- 
alyzed the ARPES data and also showed strong coupling 
features at 21.4 meV, 38.3 meV, 44.4 meV, 60.0 meV, 
and 76.2 nieV, It is remarkable that the peak positions 
at 44.4 meV, 60.0 meV, and 76.2 meV for BSCCO match 
precisely with the peak positions at 45 meV, 60 meV, and 
76 meV for underdoped La2_xSr:rCu04 (Ref. This 
consistency suggests that the fine structures seen in the 
high-resolution ARPES data are indeed related to strong 
electron-phonon coupling to multiple phonon modes. 

Strong electron-phonon coupling to multiple phonon 
modes have been also seen in several high-quality single- 
particle tunneling spectra [Hol - fssj . For conventional 
superconductors, strong electron-phonon coupling fea- 
tures cle arly show up in single-particle tunneling spec- 
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32, 54 



55|. The energies of the phonon modes 
coupled to electrons can be precisely determined from 
tunneling spectra. More specifically, the energy posi- 
tions of the peaks in —d^I/dV^, measured from the 
isotropic s-wave superconducting gap A, correspond to 
those of the peaks in the electron-phonon spectral fnnc- 
tion a^{Lo)F{Lo) (Refs. [H, 0, HI] ) . Therefore, if the cou- 
pling strength a^(cj) changes smoothly with uj, the fine 
structures in the phonon density of states F{uj) will have 
a one-to-one correspondence to those in —d^I/dV^. Such 
a conventional approach to identify the electron-boson 
coupling features would have been extensively applied 
to high- Tc superconductors if the superconducting gap 
were isotropic. Since the superconducting gap is highly 
anisotropic in high- Tc superconductors, it is difficult to 
reliably determine the energies of bosonic modes if tun- 
neling current is not directional. This may explain why 
the electron-boson coupling structures extracted from 
earlier tunneling spectra are not reproducible among dif- 
ferent groups [Ho, (ill, [5^. 

Zhao [521 analyzed high-quality tunneling spectra of 
YBa2Cu307_5 and Bi2Sr2CaCu208+5. The spectra of 
the second derivative of tunneling current d^//dT^^ in 
both compounds show clear dip and peak features due 
to strong coupling to the bosonic modes mediating elec- 
tron pairing. The energy positions of nearly all the peaks 
in — d^//dy^-like spectra match precisely with those in 
the phonon density of states obtained by inelastic neu- 
tron scattering. Similar results have also been reported 
in a Lai.84Sro.i6Cu04 thin film [53]. These results con- 
sistently demonstrate that the bosonic modes mediating 
the electron pairing are phonons and that the large at- 
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tractive interaction should arise primarily from strong 
coupling to multiple phonon modes. 
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FIG. 2: The -d^(ReE)/dcj^ spectrum (S is electron self- 
energy) and the —Ag/Au) tunneling spectrum of BSCCO crys- 
tal. After Ref. 



In Fig. 2, we compare the strong coupling features re- 
vealed in the tunneling spectrum and in the electron self- 
energy spectrum of BSCCO. The spectrum of the second 
derivative of the real part of electron self-energy (S) is 
reproduced from Ref. [49[ . It is striking that the peak fea- 
tures in (— d^(ReI])/da;^) match very well with those in 
(— A(d5/dcLi), where g is the renormalized tunneling con- 
ductance proportional to the first derivative of tunneling 
current [52|. This excellent agreement between the tun- 
neling and self-energy spectra further demonstrates that 
the observed strong-coupling features in both spectra are 
intrinsic. 

On the other hand, 
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mis- 



the authors of Ref. 
takenly assigned the energy (10.5 meV) of a peak 
position in -|-d^//dV^^ spectra of an electron-doped 
Pro.88LaCeo.i2Cu04 to the energy of a bosonic mode. 
Since this mistakenly assigned mode energy (10.5 meV) 
is very close to the energy (9.5-11 meV) of the mag- 
netic resonance-like mode measured by inelastic neutron 
scattering 58, [H^, the authors [HI] conclude that the 



magnetic resonance mode mediates electron pairing in 
electron-doped cuprates. Zhao (gqI reanalyzed the tun- 
neling data and found strong coupling to a bosonic mode 
at about 16 meV, in quantitative agreement with early 
tunneling spectra of Ndi.85Ceo.i5Cu04 (Ref. [6l[). The 
correctly assigned mode energy of about 16 meV rules out 
its connection to the magnetic resonance mode which has 
energy of 9.5 meV in NCCO (Ref. [Hi]) and 11 meV in 
PLCCO (Ref. (sij). 

The pairing mechanism based on strong coupling to a 
magnetic resonance mode predicts a pronounced dip fea- 
ture at the energy of the magnetic resonance mode in tun- 
neling conductance spectra and angle-resolved photoe- 
mission spectra [62] . Some tunneling spectra of optimally 



doped BSCCO appear to support this prediction. For ex- 
ample, the energy of the dip feature in the tunneling spec- 
trum of BSCCO is found to be about 45 meV (Ref. 
close to the magnetic resonance energy of about 43 meV 
(Ref. [oil). However, the inversion of the same tunnel- 
ing spectrum [sif yields a very strong coupling feature 
at about 20 meV, in quantitative agreement with the 
mode energies of several overdoped BSCCO crystals in- 
ferred from the electron self-energies along the antinodal 
direction [(mJ. Furthermore, this theoretical prediction 
also contradicts more recent atomic-resolution tunneling 
spectra in various BSCCO crystals. The observed dip 
features in these crystals center around 26 meV and are 
nearly independent of doping or Tc [s^j . Since the energy 
of the magnetic resonance mode is found to be propor- 
tional to Tc (Ref. [oil), the fact that the dip features 
are nearly independent of doping suggests that they are 
not caused by strong coupling to the magnetic resonance 
modes. Furthermore, the energy (26 meV) of the dip 
feature in optimally doped BSCCO is too low compared 
with the energy (43 meV) of the magnetic resonance 
mode and the energy (35-40 rneV) of the out-of-phase 
bond buckling oxygen modes . 

Further evidence for no magnetic coupling is magneto- 
optical experiments on various cuprates [68t . .69,] . The 
data have shown that the electron-boson spectral func- 
tion is independent of magnetic field [68*1 , in contradiction 
with the theoretical prediction based on the magnetic 
pairing mechanism (see Fig. 9 of Ref. 68] ) . 



CO-EXISTENCE OF BIPOLARONS (FORMED 
BETWEEN APICAL AND IN-PLANE OXYGEN 
HOLES) WITH IN-PLANE OXYGEN HOLE 
POLARONS 



It is known that undoped parent cuprates are charge- 
transfer insulators, so doped holes should mainly reside 
on the oxygen orbitals. Bulk-sensitive x-ray-absorption 
experiment on (Yi_2:Caa;)Ba2Cu307_j, confirms that 
doped holes mainly reside on both in-plane oxygen and 
apical oxygen orbitals [to'I . What is striking is that doped 
holes via oxygen doping in the CuO chains are almost 
equally distributed to the apical and in-plane oxygen or- 
bitals 1?^. The Ca substitution for Y yields hole doping 
into the in-plane oxygen orbitals only and no supercon- 
ductivity is seen in a compound with 10% in-plane oxygen 
holes and negligible apical oxygen holes. This remark- 
able experiment clearly demonstrates that the apical oxy- 
gen holes are crucial to high-temperature superconduc- 
tivity. The data provide strong support for the theory of 
bipolaronic superconductivity proposed by Alexandrov 
and Mott i^]. The theory is based on mobile intersite 
bipolarons that are formed between apical and in-plane 
oxygen holes. Different types of bipolarons in La2Cu04 
were investigated by Catlow et al. 



7ll ] with computer 



simulation techniques based on the minimization of the 
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ground-state energy without the hopping term. An in- 
tersite bipolaron formed between an in-plane and an api- 
cal oxygen-hole polaron (denoted by AP bipolaron) was 
found energetically favorable with the binding energy of 
0.14 eV. An intersite bipolaron formed between two in- 
plane oxygen hole polarons (denoted by PP bipolaron) 
is also possible but with a much smaller binding energy 
(0.06 eV). Both electron-phonon interactions and bipo- 
laron binding energies will decrease when the hoping term 
is turned on. Doping will further reduce the bipolaron 
binding energy due to the screening of charge carriers. 
Therefore, it is possible that only AP bipolarons are sta- 
ble in doped cuprates. For YBCO, about half of the 
inplane oxygen holes will be combined with the apical 
oxygen holes to form intersite bipolarons while the re- 
maining in-plane oxygen holes are mixed with in-plane 
Cu holes and may be bound into fc-space pairs below Tc. 
For LSCO, the number of the apical oxygen holes might 
be quite close to that of the in-plane oxygen holes since 
there are two apical and two in-plane oxygen atoms per 
formula unit. 

An important proof for this simple model is pro- 
vided by the data of the normal-state susceptibility x{T) 
(Ref. [tI, [t^]). The temperature dependence of the 
normal-state susceptibility in cuprates is very different 
from conventional metals. It was found 72] that x{T) 
exhibits a broad peak at a temperature Tmax- Such be- 
havior was interpreted in terms of the 2-dimensional (2D) 
antiferromagnetic correlation among Cu spins. However, 
Nakano et al. [73| have shown that, at low tempera- 
tures, xiT) strongly deviates from the prediction of the 
2D AF model. Furthermore, the magnitude of xi^) 
is proportional to the electronic specific heat [t^], im- 
plying that the xi^) reflects the density of states of 
the conduction electrons rather than the localized spin 
correlation. Alternatively, Alexandrov, Kabanov and 
Mott (AKM) explained the normal-state susceptibility 
data on the basis of their small (bi)polaron theory [75| 
which predicts a temperature-dependent spin suscepti- 
bility: xakm{T) = B.^T~2 exp(-A6p/2T), where A^p 
is the bipolaron binding energy and Boo is a constant 
depending on the effective masses of polarons and bipo- 
larons [75| . If we consider possible coexistence with in- 
plane oxygen-hole polaron Fermi-liquid, the total suscep- 
tibility for the two-component system [AP (bi)polarons 
-I- in-plane oxygen- hole polarons] is xC^) = fsXAKM{T) 
+ (1 — fs)XF + XcV- Here fs is the fraction of the AP 
(bi)polarons, xf is the susceptibility contributed from 
the Fermi-liquid-like polarons, and XcV is the total or- 
bital contribution. Then 

xiT) - fsB^T-i exp(-Abp/2fcsT) + Xo (8) 

Fig. 3 shows the temperature dependence of the 
normal-state susceptibility for Lai_a;Sra;Cu04. For x = 
0.20, a small Curie paramagnetic susceptibility has been 
subtracted. The lowest solid line marks the total orbital 
contribution, as determined from the ^^Cu Knight shift 
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FIG. 3: The temperature dependence of the normal-state sus- 
ceptibility for Lai_2,Srj,Cu04 (reproduced from Ref. _76]). 
The lowest solid line marks the total orbital contribution, 
as determined from the ^^Cu Knight shift of Lai_a;Sra;Cu04 
(Ref. Ill) 

of Lai_:,Sr:,Cu04 (Ref. [781). The solid lines in Fig. 3 
represent the fitting curves using Eq. [S] The fitting is 
excellent for all the compositions. It is interesting to see 
that for X > 0.09, (1 — fs)xF 7^ 0, that is, the Fermi- 
liquid type carriers are also present. From the more ex- 
tended x(T) data [tJ], one sees that for x < 0.10, only 
(bi)polaronic charge carriers exist, and that the fraction 
of the Fermi-liquid carriers increases monotonically with 
a; for a: > 0.10. Therefore, the theory of bipolaronic su- 
perconductivity should be applied for x < 0.10. For x > 
0.1, it is possible that, due to a decrease of the charge- 
transfer gap, more in-plane Cu holes are introduced so 
that the component of the Fermi-liquid-like carriers in- 
creases substantially. In the overdoped region (e.g., x 
= 0.2), the majority of the carriers are polarons at Tc, 
so superconductivity should arise from fc-space pairing of 
Fermi-liquid-like polarons. 

Now we discuss how this simple model is also compati- 
ble with ARPES resuhs [zi] for Lai_:rSr:rCu04. The hole 
pocket centered around (7r/2, 7r/2) has been observed by 
ARPES for X < 0.18. The hole pocket should be asso- 
ciated with the mixed states of the inplane oxygen p or- 
bitals and Cu dx^_y2 orbitals. A very flat band below the 
Fermi energy is seen near (tt, 0) at all the doping levels. 
This band has been attributed to the AP hole bipolaron 
state [80] . Above optimal doping, the hole pocket around 
(7r/2, 7r/2) appears to be connected with the hole pocket 
around (tt, 0) to form a large Fermi surface predicted 
by local density approximation. Nevertheless, the mea- 
sured bare plasma energy for x > 0.15 is only about 1.8 
eV (Ref. 81]), which is much smaller than the predicted 
value of about 2.9 eV (Ref. [s^). Therefore, the apparent 
connection between the two hole pockets may arise from 
the band-tailing phenomenon in doped semiconductors 
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THE INTRINSIC PAIRING SYMMETRY IN THE 
BULK OF CUPRATES 

An unambiguous determination of the pairing symme- 
try in cuprates is crucial to the understanding of the pair- 
ing mechanism of high-temperature superconductivity. 
Many experiments have been designed to test the pair- 
ing symmetry in the cuprate superconductors. However, 
contradictory conclusions have been drawn from differ- 
ent experimental techniques I84H90I. 92-981, which 
can be classified into being bulk-sensitive and surface- 
sensitive. For example, the magnetic penetration depth 
measurements and polarized Raman scattering experi- 
ments are bulk-sensitive. Angle-resolved photoemission 
spectroscopy is in general a surface-sensitive technique. 
However, the ARPES data for Bi2Sr2CaCu208+i, should 
nearly refiect the bulk properties since the cleaved top 
surface contains an inactive Bi-0 layer, and the su- 
perconducting coherent length along the c-axis is very 
short. The single-particle tunneling experiments along 
the Cu02 planes can probe the bulk electronic density 
of states since the mean free path is far larger than 
the thickness of the degraded surface layer [99|. In 
contrast, the phase-sensitive experiments based on the 
planar Josephson tunneling are rather surface sensitive, 
so that they might not probe the intrinsic bulk super- 
conducting state if the surfaces are strongly degraded. 
Therefore, the surface- and phase-sensitive experiments 
do not necessarily provide an acid test for the intrinsic 
bulk gap symmetry. 

Zhao [100] has identified the intrinsic bulk pairing sym- 
metry for hole-doped cuprates from the existing bulk- 
and nearly bulk-sensitive experimental results such as 
ARPES, magnetic penetration depth, single-particle tun- 
neling spectra, and nonlinear Meissner effect. These 
experimental results consistently show that the domi- 
nant pairing symmetry in hole-doped cuprates is an ex- 
tended s-wave with eight line nodes, that is A{6) — 
As -|- AgC0s4^, where the g component Ag is larger 
than the s-wave component A^, and ^ is measured from 
the Cu-0 bonding direction. Zhao jl01| has also stud- 
ied the intrinsic pairing symmetry for optimally electron- 
doped cuprates from bulk-sensitive data of Raman scat- 
tering, optical conductivity, magnetic penetration depth, 
directional point-contact tunneling spectra, and nonmag- 
netic pair-breaking effect. He has shown that all these 
bulk-sensitive data are in quantitative agreement with a 
nearly isotropic s-wave gap. The bulk-sensitive specific 
heat data of optimally electron-doped cuprates are also 
in quantitative agreement with a nearly isotropic s-wave 
gap [l02. 

In particular, bulk and phase-sensitive nonmagnetic 
pair-breaking effects can make clear distinction a mong 
different gap symmetries. Numerical calculations 101 
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FIG. 4; Tc as a function of residual resistivity in optimally 
doped Ndi,g5Ceo.i5Cu04_y. The solid line is numerically cal- 
culated curve in terms of any d-wave gap and the dotted line is 
the calculated curve for an s-wave gap function proportional 
to (1 - 0.11 cos 461). After Ref. [lOll ]. 



for impurity pair-breaking effects indicate that Tc can 
be suppressed to zero very rapidly with both magnetic 
and nonmagnetic impurities for a c?-wave gap: A{9) = 
Ad cos 20 or a g-wave gap. In contrast, Tc is hardly sup- 
pressed wit h non magnetic impurities for nearly isotropic 
s-wave gap 10ll | . The fact that Tc is nearly independent 
of the residual resistivity in optimally electron-doped 
Ndi.85Ceo.i5Cu04_y (see Fig. 4) should r ule g ut any d- 
wave gap including a nodeless c?-wave gap |l03l | . 
Recent 



calculation based on t — J model [104 1 
shows that the nonmagnetic pair-breaking effect be- 
comes less effective due to strong electron-electron cor- 
relation. This would explain the insensitivity of Tc 
to the residual resistivity in the optimally electron- 
doped Ndi.85Ceo.i5Cu04_y . However, this cannot con- 
sistently explain the very rapid Tc suppression with 
the residual resistivity in a hole-doped Lai.8oSro.2oCu04 
(Fig. 5b) unless the electron-electron correlation in 
Lai.8oSro.2oCu04 is negligibly small compared with 
that in Ndi.85Ceo.i5Cu04_j,. Since the renormal- 
ized plasma energy HQ* = 1.00 eV (see below) of 
Lai.8 oSrn 20CUO4 is even much lower than that (1.64 



eV) jlOSl ] for Ndi.85Ceo.i5Cu04_y, the correlation ef 



feet of the former should be much stronger than that 
of the latter. If the t — J model were relevant, Tc of 
Ndi.85Ceo.i5Cu04_y would be suppressed more rapidly 
than that of Lai.8oSro.2oCu04. This is in total contra- 
diction with the experimental results (see Figs. 4 and 5). 

In fact, we can quantitatively explain the magnetic 
penetration depth and non-magnetic pair breaking effect 
in Lai.8oSro.2oCu04 even without taking into account 
electron-electron correlation. Fig. 5a shows the temper- 
ature dependence of the superfluid density [c>c 1/A^^(T)] 
for Lai.8oSro.2oCu04. The solid line is the numerically 
calculated curve for an extended s-wave (s*-wave) gap 
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FIG. 5: a) Temperature dependence of the superfluid density 
[oc l/X^b{T)] for Lai.goSro.2oCu04. The solid line is the nu- 
merically calculated curve for an extended s-wave gap func- 
tion: A(6') = 6.38(0.24 -|- cos 40) meV. The data are taken 
from Ref. [l06l ]. b) Tc dependence on the residual resistivity 
for Lai.goSro.2oCu04. The solid line is calculated curve for 
the extended s-wave gap function: A(6) = 6.38(0.24-1- cos 40) 
meV. The data are taken from Ref. IIOTH. 



function: A{9) = 6.38(0.24 -|- cos 46*) meV using the 
standard formula |lOCl| and the BCS temperature de- 
pendence of the gap. This gap function has a simi- 
lar form as that for YBa2Cu307 with a similar dop- 
ing |lOOl |. Fig. 5b shows Tc dependence on the resid- 
ual resistivity for Lai.8oSro.2oCu04. The solid line is 
calculated curve for the extended s-wave gap function: 
A{e) = 6.38(0.24 + cos 46*) meV obtained from the pene- 
tration depth data (Fig. 5a) using the standard formula 
[lO0| | and the renormalized plasma energy hVl* = 1.00 eV 
that is calculated from the measured zero-temperature 
penetration depth (2000 A) (see Fig. 5a). It is apparent 
that the calculated curve is in quantitative agreement 
with the data. In contrast, the calculated curve (dashed 
line) for d-wave gap is significantly off from the data. 

When hole doping is reduced, the s-wave component 
increases so that the nonmagnetic pair-breaking effect 
becomes weaker. Since the bipolaronic charge carriers 



become dominant in underdoped cuprates, nonmagnetic 
impurities do not break the pairs, but can still suppress 
Tc via localization effect. 

With the two-carrier scenario in hole-doped cuprates, 
we can consistently explain surface and phase-sensitive 
experiments which probe d-wave order-parameter sym- 
metry at surfac es and int erfaces. Experiments on hole- 
doped cuprates (l08l Il09| have shown that surfaces and 
interfaces are significantly underdoped so that bipola- 
ronic charge carriers are dominant over polaronic Fermi- 
liquid-like ones. Since the order-parameter symmetry 
of B ose-Einstein condensate of bipolarons are d-wave 
[llOj l , the surface and phase-sensitive experiments should 
probe the dominant d-wave ord er-p arameter symmetry, 
in agreement with experiments [92| . For electron-d oped 
(n-type) cuprates, the penetratioi i depth l llll . Ill2| and 
point-contact tunneling spectra |ll3l4ll7f consistentlv 
suggest that the gap symmetry in deeply underdoped 
samples should be d-wave and change to a nodeless s- 
wave when the doping level is above a critical value, 
in agreement with the theoretical prediction based on a 
phonon- mediated pairing mechanism [118]. This scenario 
can also explain the d-wave gap symmetry inferred from 
surface-sensitive experiments if surfaces or interfaces are 
deeply underdoped. Another possibility is that supercon- 
ductivity in deeply underdoped rt-type cuprates is also 
due to the Bose-Einstein condensation of local pairs. In 
this case, the superconducting qua si-particle gap is found 
to have a d-wave symmetry [119j . In sharp contrast to 
the standard theories of Bogoliubov quasi-particle exci- 
tations, the quasi-particle gap is shown [119] to originate 
from anomalous kinetic process, completely unrelated to 
the pairing symmetry. 

CONCLUDING REMARKS 

We have presented some crucial experiments that place 
essential constraints on the pairing mechanism of high- 
temperature superconductivity. The observed uncon- 
ventional oxygen-isotope effects in cuprates have clearly 
shown strong electron-phonon interactions and the exis- 
tence of polarons and/or bipolarons. Angle-resolved pho- 
toemission and tunneling spectra have provided direct 
evidence for strong coupling to multiple-phonon modes. 
In contrast, these spectra do not show strong coupling 
features expected for magnetic resonance modes. Angle- 
resolved photoemission spectra and the oxygen-isotope 
effect on the antiferromagnetic exchange energy J in un- 
doped parent compounds consistently show that the po- 
laron binding energy is about 2 eV, which is over one 
order of magnitude larger than J = 0.14 eV. In addi- 
tion, a photoinduced latt ice e xpansion experiment |l20| 
and various optical data [l2l'| also provide evidence for 
strong electron-phonon coupling. 

The normal-state spin-susceptibility data of hole- 
doped cuprates indicate that intersite bipolarons are the 
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dominant charge carriers in the underdoped region while 
the component of Fermi-hquid-hke polarons is dominant 
in the overdoped region. This two component scenario is 
also supported by the femtosecond optical spectroscopic 



experiment on (Yi_j;Ca2:)Ba2Cu30 
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Various bulk- 



sensitive experiments consistently demonstrate that the 
intrinsic gap (pairing) symmetry for the Fermi-liquid- 
like component is anisotropic s-wave while surface- and 
phase-sensitive experiments see d-wave order-parameter 
symmetry at intrinsically underdoped surfaces and inter- 
faces where the dominant bipolaronic charge c arriers ar e 
condensed into a d-wave superconducting state lid 119{ . 



gzhao2@calstatela.edu 
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